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Abstract 

The physical representation of the general double-Reissner-Nordstrom solution is obtained by 
rewriting the N = 2 Breton-Manko-Aguilar electrostatic solution in the Varzugin-Chistyakov 
parametrization (Mj, Qi, R). A concise analytical formula is derived for the interaction force 
between two arbitrary Reissner-Nordstrom constituents, and an example of the equilibrium con- 
figuration involving two oppositely charged particles which confirms earlier Bonnor's prediction of 
the existence of such configurations is given. 



PACS numbers: 04.20.Jb, 04.70.Bw, 97.60.Lf 



I. INTRODUCTION 



The history of exact solutions of the Einstein-Maxwell equations for two aligned charged 
masses of the Reiss„er-Nordstr5m type 0. 3 begins in 1917 wth the two-body e.eetro- 
static solution of Weyl's class [3[ which is able to describe constituents whose masses Mj 
and charges Qj, z = 1,2, satisfy the relation M1Q2 — M2Q1 = 0. When the parameters Qi 
are defined by the equalities Qi = ±Mi, Q2 = ±M2,^he resulting solution belongs to the 
Majumdar-Papapetrou family of extreme black holes P, ^, and in that case the constituents 
are in equilibrium which is independent of the coordinate distance R separating them. The 
study of the general exact model characterized by arbitrary Mj and Qi was started only a 
decade ago by Perry and Cooperstock who adjusted a known stationary axisymmetric 
electrovac solution [7] for investigating the double-Reissner-Nordstrom equilibrium prob- 
lem. The solution analyzed in 6|] was not presented in a closed analytical form since the 
parameters a„ it contained were implicit functions of the parameters of the axis data. A year 
later the multi-soliton electrostatic solution |^ representing a system of aligned Reissner- 
Nordstrom particles was constructed with the aid of SibgatuUin's integral equation method 
91] in a concise analytical form thanks to the use of the objects a„ as arbitrary parameters. 
Its = 2 specialization, henceforth referred to as the BMA solution, fully describes the 
general 5-parameter double-Reissner-Nordstrom spacetime, the corresponding subclass of 
equilibrium configurations being defined by a very simple balance equation. 

In a recent paper [lo| Alekseev and Belinski have been able to parametrize the 4- 
parameter subfamily of the BMA solution representing two Reissner-Nordstrom particles 
in equilibrium in terms of the Komar masses and charges of the constituents, the coordinate 
distance separating the balancing particles being a function of the mass and charge param- 
eters. Since no technical details of the derivation have been provided in [l^, it is one of 
the motivations of the present paper to demonstrate that Alekseev and Belinski's results are 
obtainable straightforwardly, and by purely algebraic manipulations, as a particular case of 
the formulas of the paper 8| rewritten in the parametrization discovered and elaborated for 
the double-Reissner-Nordstrom problem by Varzugin and Chistyakov HI]. For instance, for 
obtaining Eq. (10) of [l^, one only needs to rewrite the coefficient 012 of the BMA solution 
(see Eqs. (2) and (9) below) in the new parameters. 

However, the removal of mystery from Alekseev and Belinski's article is not the main 



objective of our research. The present paper has as its principal goal the presentation of 
the entire 5-parameter BMA family of electrostatic solutions in the Varzugin-Chistyakov 
parameterization and derivation on its basis of the formula for the interaction force between 
two arbitrary Reissner-Nordstrom constituents. The extreme technical complexity of this 
task probably explains why neither Varzugin and Chistyakov themselves nor later on Alek- 
seev and Belinski have been able to accomplish it. Fortunately, the difficulties of the analytic 
computer processing which at first glance look insuperable can be circumvented with the aid 
of some special tricks, but to give an idea about the scale of these difficulties, it is sufficient 
to mention that for instance an attempt to directly rewrite and then factorize the coefficient 
ai3 in terms of the new parameter set exhausts the memory of a 12 Gb RAM computer. 

The paper is organized as follows. In the next section the BMA electrostatic solution in 
its original 'a — /5' parametrization is briefly reviewed. Sec. Ill is devoted to the Varzugin- 
Chistyakov parametrization and its relation to the parameters of the BMA solution. In 
Sec. IV the general double-Reissner-Nordstrom spacetime is written in terms of the coor- 
dinate distance R and physical parameters Mj and Qi, after which the derivation of the 
desired formula for the interaction force is carried out. Concluding remarks are given in 
Sec. V. 



II. THE BMA SOLUTION IN THE 'a - /?' PARAMETRIZATION 



The Ernst potentials 8 and $ [12| of the BMA solution describing two aligned Reissner- 
Nordstrom particles, and the corresponding metric functions f{p,z), ^[p, z) entering the 
static axisymmetric line element 

ds^ = f-'[e^^{dp^ + dz^) + pMif^] - fdt^ (1) 

where p and z are Weyl's cylindrical coordinates, are defined by the formulas 
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In the above formulas the arbitrary parameters are q;„, n — 1,4, and / = 1,2, which 
can take on arbitrary real values or occur in complex conjugate pairs, all other constant 
quantities, i.e., D(q;„), and /i;(q;„), being defined in terms of and /3; through the 

relations 



/i 



+ 



h 



D{an) = {an- Pi){an- P2), f{an)^ r, ■ r, . 

an - Pi an- P2 

hi{an) = ei + 2 f if {an), h2{an) = 62 + 2/2/(q;„), 



(3) 



where 



ei = 



62 



[-2S4 + {Pi + P2)S3 - 2/3i/?252 + /?i'(3/?2 - Pl)si 

+2(3!{(3i - 2(32)]{(3i - (32)-' - 2/i/2(/3i - (32)-\ 

[2s4 - {(3i + (32)ss + 2(3i(32S2 - Pi{SPi - (32)si 

-2Pl{(32 - 2A)](/3i - P2)-' + 2/i/2(/9i - P2)-\ 
4 

^Qij, S2 ^ ttittj, 83-.= ^ aiajak, S4 :^ aia2a3a4^, (4) 

l<i<j<4 l<i<j<k<4 



i=l 



and 



fi 



/I 



nn=l(/^2 - an) 



(A-W ' {Pi -(32)^ ■ 

On the upper part of the symmetry axis the potentials S and $ behave themselves as 
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and it is worth mentioning that Ci, 62, /i, /2, together with the simple poles Pi and P2, 
can in principle be used as arbitrary parameters of the solution, but then the parameters 
«„ would be implicit functions of e;, fi, Pi defined as roots of an algebraic quartic equation 
(see js] for details), and in such a case one would need to estimate them numerically in 
applications. 

In view of the invariance of the field equations with respect to an arbitrary constant shift 
Zq along the 2;-axis, the set {an, Pi} consisting of 6 parameters can define only 5 physical 
characteristics of the binary system since the constants a„ determining the location of sources 
on the symmetry axis (see Fig. 1) can be always set to satisfy, say, the equation ^ a„ = 0. 

The total mass M and total charge Q of the system are given by the expressions 

M = -i(ei + e2), Q = U + h, (7) 

and these formulas play an important role in the reparametrization which will be carried 
out in the next section. 

Although the BMA solution is analytically extended and therefore describes any combi- 
nation of the black-hole and hyperextreme constituents, in what follows we shall be mainly 
concerned with the case of two black holes because of its greater physical importance, making 
where necessary only brief explanatory remarks about the binary systems involving hyper- 
extreme constituents. In the subextreme case, all a„ are real constants to which one can 
assign, without loss of generality, the order 

«! > ^2 > «3 > «4, (8) 

and the locations of the upper and lower constituents are defined, respectively by the seg- 
ments p = 0, 0^2 < -2 < ai and p = 0, 0:4 < z < 0:3 of the symmetry axis (see Fig. 1). The 
part p = 0, as < z < ^2 of the axis represents a Weyl strut Q which prevents the 
constituents from falling onto each other. The strut is absent in the cases when the metric 
function 7 is zero on the above mentioned interval separating the constituents and, as was 
shown in 8|, this happens when the parameters of the BMA solution satisfy the equation 
ai2 + '334 = 0. The latter equation, as was already remarked in [l^, can be further trivialized 
by observing that 034 = ai2, so that the condition of equilibrium of two Reissner-Nordstrom 
constituents due to the balance of the gravitational and electrostatic forces is defined by 
vanishing of the coefficient ai2 alone: 

"12 = <^ {a'i- p2){ai- Pi)hi{a'i)h2{ai)-{a^- Pi){ai- p2)hi{ai)h2{a'i) =Q. (9) 



By assigning particular values to ai, a^-, as, 04, /5i and finding numerically the corre- 
sponding value of /52 from Eq. ([9]), the equilibrium configurations between a sub- and a 
hyperextreme constituent were found in [sj], thus confirming Perry and Cooperstock's earlier 



results {g] obtained with the aid of an exact solution, and Bonnor's analysis of the approxi- 
mate charged two-body problem 16|. Although the Komar masses and charges of the 
balancing constituents were also calculated in [s^], it was not observed there that the 
parametrization is ideally appropriate for the introduction of the Komar quantities explicitly 
into the double-Reissner-Nordstrom spacetime. 



III. THE VARZUGIN CHISTYAKOV PARAMETRIZATION 

I — I 

Following the main ideas of Varzugin's study of the stationary vacuum case p^] , in the 



paper 



11| Varzugin and Chistyakov extended Carter's analysis of a single black hole pi 
to the system of aligned charged rotating black holes supported by struts. They suc- 
ceeded in relating their boundary-value problem, parametrized in terms of the individual 
Komar quantities and some additional parameters such as angular velocities of the horizons 
or electric potentials, to the matrix Riemann-Hilbert problem which they were solving on 
the z-axis. They used the properties of the reduced monodromy matrix T{}z) to obtain 
the non-linear algebraic constraints on the parameters of the boundary- value problem, thus 
introducing the independent parameters in terms of which all other parameters might be 
expressible. Varzugin and Chistyakov applied their approach to the case of two arbitrary 
Reissner-Nordstrom black holes and obtained the expressions for black holes' irreducible 
masses in terms of the individual Komar masses Mj, Komar charges Qi and the coordi- 
nate distance R between the constituents (5 independent parameters in total), formulas 
(46) constituting the central result of the paper As will be seen below, the knowl- 

edge of irreducible masses in terms of Mj, Qi and R is sufficient for reparametrizing the 
BMA solution in terms of the latter quantities, so it will be correct to call the parameter 
set {Ml, M2, Qi, Q2, -R} the Varzugin-Chistyakov parametrization of the double-Reissner- 
Nordstrom problem, even though Varzugin and Chistyakov have never taken an interest in 
solving that problem outside the symmetry axis. 

The Varzugin-Chistyakov formulas for the irreducible masses Oi written in Alekseev and 



Belinski's manner [lO|] have the form 11 1 



and we refer the reader to the paper ll|] for all the details of their derivation. Note that the 
irreducible masses are denoted here as cTj instead of Varzugin and Chistyakov's because 
these quantities take pure imaginary values in the hyperextreme case and hence should not 
be confused with the genuine masses. 

Formulas ffTOl) provide us with the immediate reparametrization of the constants a„ of 
the BMA solution: 

ai = jR + a2, 0:2 = |i? — (X2, as = —^R + ai, 04 = — |-R — (Xi, (11) 

because, by definition, 

(Ti = |(q;3 — 04), cT2 = |(ai — 02), -R = + ^2 — as — «4), (12) 

where R is the coordinate distance between the centers of the black hole horizons (see 
Fig. 1). Mention that we have chosen the origin of the coordinate system in such a way 
that V an = 0, and we place the constituent denoted by index 2 above the constituent with 

nn 

index 1 to meet the conventions of papers [10, lUj • 

Our next objective is to express the remaining parameters of the BMA solution, the poles 
/3i and (32, in terms of Mj, Qi and R. This can be done in the following way. First, the 
substitution of ei and 62 from (jl]) into the first formula in ([7]) yields the equation 

A + /?2 = |si - M = -M, (13) 

because si = by virtue of our choice of a„. Then, denoting the right-hand sides of the 
first and second equalities in ([5]) as uji and lj2, respectively, we pass from the second formula 
in ([7]), namely, Q = /i + /2, to the equation 

{Ui +UJ2- Q'^f - 4:LUiUJ2 = 0, (14) 

the left-hand side of which, after substituting into it the explicit expressions for a„, Ui, U2, 
and also (32 from (|T3i) . factorizes into a pair of polynomials quadratic in (3i. The solution 
which is consistent with the second equality in ([71) is 



/3i = -|(M + D), D = y/R^ + (Ml - M2)(Mi - M2 + 2R) + AQ2{Qi - 2/x), (15) 



where M = Mi + M2. The expression for (32 is now readily obtainable from fll3p . yielding 

/?2 = -i(M-D). (16) 

Once the constants a„ and f3i are reparametrized in terms of Mj, Qi and R, no difficulty 
arises in rewriting the quantities /; and ei defined by formulas and (jlj) in terms of the 
new parameter set. Observe that from ([5]) and ([7]) follows that 

/l -/2 =Q(/l-/2)=^l-^2 =^ h- f2 = {uJi-U2)Q-\ (17) 

SO that the expressions for fi are most easily obtainable through the formulas 

h = \[Q + {ui - u2)Q-\ /2 = |[Q - (^1 - oj2)Q-\ (18) 

giving as the result 

h = [QD + (Ml - M2 + R){Qi - Q2) + AMiQ2]{2D)-\ 

f2 = [QD - (Ml -M2 + R){Qi - Q2) - miQ2]{2D)-\ (19) 

where Q = Qi + Q2- 

For ei we readily get from (jlj): 

ei = -[M(M + D) + R{Mi - M2) - 2fiQ]D-\ 

62 = [M(M -D) + R{Mi - M2) - 2fiQ]D-\ (20) 

Therefore, we have obtained all the necessary formulas for being able to rewrite the BMA 
solution in the Varzugin-Chistyakov parametrization. 

IV. THE BMA SOLUTION IN THE PHYSICAL PARAMETRIZATION AND 
THE FORMULA FOR THE INTERACTION FORCE 

The simplest application of the formulas obtained in the previous section is rewriting the 
BMA equilibrium condition (j9j) in terms of the new parameters. The substitution of (j3j), 
(jTOj), (jll]), (jlS]), (jlSj) and (j20j) in Eq. © immediately leads to 



MiM2-(gi-/i)(g2 + /i) = o, 



(21) 



and no computational problem arises during this calculation. One easily recognizes in fl2ip 
the balance equation (10) of the paper flo| . 

Turning now to the presentation of the reparametrized general BMA solution, it can 
be remarked that the main technical difficulty during the rewriting of the coefficients bi 
and Ci from ([2]) was finding their concise expressions which would replace the cumbersome 
intermediate formulas. The coefficients ai2 and 034 which enter the expression for A are 
the simplest ones for working out (besides, 034 = a^) and are obtainable at once in their 
final form. On the other hand, we did not succeed in a straightforward evaluation of the 
coefficients ai3, 014, 023 and 024; the only reasonable way of finding their reparametrized form 
was to process them by smaller fractions. The reparametrized BMA solution eventually takes 
a concise form by introducing the constant objects u and n defined by 

u:=R^ -al-al + 2fi\ k := M,M2 - {Qi - fi){Q2 + f^), (22) 

with which formulas ([2]) of Sec. II take the following elegant final form in terms of the 

physical parameters Mj, Qi and R: 

A-B _ C _ A^-B^ + C 27 _ - + C 

A + B' ~ A + B' ^~ {A + By ' ^ ~ 16afai{u + 2K,yrir2rsn 
A = 0-10-2 [z^(n + r2){r3 + r4) + AK{rir2 + r3r4)] - (/iV - 2K^)(ri - r2)(r3 - r4), 

B = 2(Tia2[(z/Mi + 2/tM2)(ri + r2) + {uM2 + 2/tMi)(r3 + n)] 

- 2ai[uii{Q2 + /i) + 2k{RM2 + liQi - li'^)]{ri - r2) 

- 2a2[iyfiiQi - fi) - 2k{RMi - 12Q2 - fi^)]irs - u), 

C = 2aicx2{HQi -fi) + 2k{Q2 + /i)](ri + r2) + \v{Q2 + /x) + 2k{Q^ - /i)](r3 + r^)} 

- 2(Ti[/iz/M2 + 2/€(/iMi + RQ2 + yui?)](ri - rs) 

- 2a2[/iz/Mi + 2/t(/iM2 - RQx + [iRW^ - r^), (23) 

where Oi are determined by formulas (fTOl) . while the reparametrized form of the functions 
r„ is 

n = ^jp'^{z-\R-a2)\ r2 = ^Jp^ + {z-lR + a2y, 

rs = ^p^ + {z + \R-ai)\ = ^ + [z + \R + a^f, (24) 

ri and r2 referring to the upper constituent endowed with Komar quantities M2, Q2, and 
and r4 referring to the lower constituent endowed with Mi, Qi. Formulas fl^2]) - fl2^ and 
(ITU]) fully describe the BMA solution in the Varzugin-Chistyakov parametrization. 



One can see that in the case of balancing constituents k = and the expressions 
simphfy further. Thus, the Ernst potentials S and $ take the form 

S = E^/E+, $ = 

E± = (Tia2(ri + r2±2M2)(r3 + r4±2Mi) 

- Hn - rs) ± 2a2{Qi - /i)][yu(r3 - n) ± 2(71 (Qa + /i)], 
C = 2aia2[{Qi - fi){ri + rg) + (Q2 + /i)(r3 + n)] 

- 2/i[M2(Ti(ri - ra) + Mia2(r3 - r4)], (25) 

while for the corresponding metric functions / and 7 we have 



i?^' 16crf cr|rir2r3r4 ' 

X [al{r3 + r,f + 4cr2(/i2 - al) - /^^(ri - r^f], (26) 

and the expressions for $, / and /~^exp(27) written in specific bipolar coordinates have 
already been given in 1^. It should be emphasized that the Ernst potentials ( |25l) satisfy 
the field equations only when the balance condition (12T1) holds, i.e., when the distance R is 
defined by the formula 

M2Q1 - MiQ2 



M and Q denoting, respectively, the total mass and total charge of the system. In view of 
this, the constant object n entering Eqs. ( l25l) and ( l26l) is given by the formula 

2(M,M2-g,g2) 

Qi - Q2 ± - 4M1M2 
and the choice of the appropriate sign in the denominator of /i depends on the analysis of a 
concrete equilibrium position. Thus, by choosing '-' in the formula (J271) . it is easy to confirm 
Bonnor's prediction about the possibility of the equilibrium of oppositely charged particles 



he made with the aid of an approximation method [16|. An example of such equilibrium 
state is the following: Mi ~ 0.2017, Qi ~ -0.2852, M2 = 1, Q2 = 2, ~ 0.4688 (the 
approximate numerical values are given up to four decimal places). The corresponding 
values for o"i and cr| are o"i ~ 0.4409, cr| ~ —1.3513, which means that the upper and lower 
constituents are, respectively, the hyperextreme and subextreme ones; besides, i? > ai, so 
that the constituents do not overlap. 



We now turn to the physically most important result which can be obtained with the aid 
of the general formulas fl23l) . As was already mentioned in Sec. II, two arbitrary Reissner- 
Nordstrom constituents are held apart by a strut when the parameters of the BMA solution 
do not satisfy the balance condition (Q. The analysis of the energy-momentum tensor asso- 
ciated with this strut permits one to introduce the interaction force between the constituents 



via the formula 



I2QI 

jr= l(e-7o_i)^ (29) 

where 70 is the constant value of the metric function 7 on the strut. In our case, formulas 
( !29|) and (!23|) conveniently give us the following expression for JF: 

^ ^ ^ ^ MiM2-{Qi-ii){Q2 + fi) . . 

y-2K R2-{M, + M2y + {Qi + Q2y' ^ ^ 

and it is surprising how remarkably simple is this result which is applicable to any pair of 
the Reissner-Nordstrom constituents. When Qi = Q2 = (the pure vacuum limit), one 
obtains fror n f[30ll the known expression for the interaction force between two Schwarzschild 
black holes [20|]: 

_ M1M2 



At large separation distances, formula fl30|) gives the Newtonian expression for the force 
between two charged particles. It should be also pointed out that the equilibrium condition 
fl2T|) is a direct consequence of the formula fl30|) if one demands JF = 0. 



V. CONCLUSION 



In the present paper we have succeeded in working out the physical representation of 
the general double-Reissner-Norstrom solution by rewriting the 5-parameter BMA electro- 
static metric in the Varzugin-Chistyakov parametrization. Formulas fl22p - (l2^ . fllOj) contain 
all known exact solutions for two non-extreme Reissner-Nordstrom constituents and are 
very suitable for the use in concrete applications. The 4-parameter subclass of the BMA so- 
lution representing the balancing constituents is a special case picked out by the equilibrium 
condition fl2Tl) . and it is worth mentioning that this subclass does not contain equatorially 
symmetric solutions because, apart from the Majumdar-Papapetrou extreme case, equilib- 
rium is only possible between a black hole and a hyperextreme object. The formula for the 



interaction force obtained in this paper is applicable to any pair of the Reissner-Nordstrom 
constituents, which is a reflection of the analytically extended character of the BMA solution. 

Since Varzugin and Chistyakov's approach to the parametrization of the multi-black- 
hole configurations has proved to be efficient in the electrostatic case, it would be also likely 
to use it, probably with some amendments, in application to the stationary axisymmetric 
electrovac solutions. As is known, the general system of aligned Kerr-Newman black 
holes is described by the Ruiz-Manko-Martm multi-soliton electrovac metric 2l|, and we 
expect that at least some of its particular cases representing the two-body configurations 
can be worked out in the physical parametrization by establishing the relationship between 
the canonical parameters of the multisoliton solution and the Komar quantities introduced 



in the paper 



111 ] via the boundary Riemann-Hilbert problem. 
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FIG. 1: The location of sources on the symmetry axis (the subextreme case). 



